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Abstract. When a rotating neutron star loses angular momentum, the reduction in the centrifugal 
force makes it contract. This perturbs each fluid element, raising the local pressure and originating 
deviations from beta equilibrium that enhance the neutrino emissivity and produce thermal energy. 
This mechanism is named rotochemical heating and has previously been studied for neutron stars 
of non-superfluid matter, finding that they reach a quasi-steady state in which the rate that the spin- 
down modifies the equilibrium concentrations is the same to that of the neutrino reactions restoring 
the equilibrium. On the other hand, the neutron star interior is believed to contain superfluid 
nucleons, which affect the thermal evolution of the star by suppressing the neutrino reactions and 
the specific heat, and opening new Cooper pairing reactions. 

In this work we describe the thermal effects of Cooper pairing with spatially uniform energy 
gaps of neutrons and protons on rotochemical heating in millisecond pulsars (MSPs) when only 
modified Urea reactions are allowed. We find that the chemical imbalances grow up to a value close 
to the energy gaps, which is higher than the one of the nonsuperfluid case. Therefore, the surface 
temperatures predicted with Cooper pairing are higher and explain the recent measurement of MSP 
J0437-4715. 
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INTRODUCTION 

The main motivation to study the thermal evolution of neutron stars is that contrasting 
theoretical predictions with the thermal emission measured from neutron stars (NSs) 
has the potential to provide constraints on their inner structure. In the existing literature, 
several detailed cooling calculations have been compared to the few estimates available 
for the surface temperatures of neutron stars (see [10] for a review and references). These 
calculations are based on the early passive cooling, which is at first neutrino-dominated. 
On the contrary, we focused our study on the late thermal evolution, where the cooling 
is driven by photon emission at ages greater than ~ lO^yr. 



Rotochemical heating 

Several mechanisms can keep NSs hot beyond the standard cooling timescale ~ lO^yr, 
among them rotochemical heating. The latter was first proposed in [4] and then improved 
in [1] by considering the internal structure of non-superfluid neutron stars via realistic 
equations of state (EOSs) in the framework of general relativity. It works as follows. 



The reduction of the centrifugal force makes the NS contract. This perturbs each fluid 
element, raising the local pressure and originating deviations from beta equilibrium, 
which are quantified by the chemical imbalances rjnpi = pin — pip — lii, where n, p, are 
neutrons and protons respectively, and / stands for leptons (electrons and muons). On 
the other hand, the neutrino reactions tend to restore the beta equilibrium, being more 
efficient as rf^pi grows. In this sense, the evolution equations for the chemical imbalances 
have the following form: 



where the terms Z„p, Z„pe, Z„pp, Wnpe, and Wnpn are constants that depend on the 
stellar structure, and is the product of the angular velocity and its time derivative 
(proportional to the spin-down power). Additionally, we have introduced the net reaction 
rate integrated over the core, defined as Ar^pi — ^n^pi — ^ pi^n- 

The evolution of the temperature of the isothermal interior, redshifted to a distant 
observer. Too, is given by the thermal balance equation [6] 



where C is the total heat capacity of the star, Ly is the total power emitted as neutrinos 
due to Urea reactions, and is the power released as thermal photons. The heating term 

produced by each Urca-type reaction is defined as = T]^gAr„pe + ri^p^M^npn- 
The most relevant feature of these equations is that, eventually, the system reaches 
a quasi-steady state where the rate at which spin-down modifies the equilibrium con- 
centrations is the same at which neutrino reactions restore the equilibrium (see figure 
1). This implies a conversion of rotational energy into thermal energy and an enhanced 
neutrino emission originated by a departure from the beta equilibrium. Thus, this mech- 
anism can keep old millisecond pulsars (MSPs) warm, at temperatures ~ 10^ K. 



Cooling curves usually consider the effects of nucleon superfluidity on the thermal 
evolution of NSs. Superfluidity is produced by Cooper pairing of baryons due to the 
attractive component of their strong interaction, and it is present only when the temper- 
ature T of the matter falls below a critical temperature T^. However, the physics of these 
interactions is rather uncertain and very model-dependent, and so is the critical temper- 
ature obtained from theory (see [3]). An important microscopic effect is that the onset 
of superfluidity leads to the appearance of a gap A in the spectrum of excitations around 
the Fermi surface. This gap in the spectrum considerably reduces the neutrino reactions 
and the specific heat involving superfluid species (neutrons and protons in the core) [9], 
and therefore, changes the evolution of rotochemical heating. 
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FIGURE 1. Evolution of the internal temperature 7U, the surface temperature Ts^oo and the chemical 
imbalances rj^pe' ^npu for ^ star with the parameters fixed to the millisecond pulsar J0437-4715, i.e. a 
mass of 1.76Mq [7], built with the A18 + dv + UIX* EOS, and a magnetic field = 2.8 • 10^ G. The 
initial conditions are T^ = \0^ K, null chemical imbalances, and an initial period of Pq = 1 ms. The error 
bar is the 90% confidence level for the surface temperature measured for the millisecond pulsar J0437- 
4715 [2] at its current spin-down parameters. Upper panel: nonsuperfluid case (null energy gaps). Lower 
panel: superfluidity of neutrons with A„ = 0.1 MeV (dashed line). 



NEUTRINO REACTIONS AND COOPER PAIRING 

We consider models in which modified Urea reactions are the main neutrino emission 
mechanism: 

n + Ni^p + Nf + e' + Ve (3) 



p + Ni + e -^n + Nf + Ve, 



(4) 



where the subscripts / and / stand for the initial and final state of the spectator nucleon 
N. 

As in [8], we compute the net reaction rate ATnpi and the neutrino emissivity for 
these reactions numerically, considering the chemical imbalances 77„p/, the energy gaps 
A in the energy spectra of the nucleons, and the temperature as free parameters (see 
details in [5]). We find that at low temperatures, i.e. T <^ Tf^pi, which is the regime of 
interest in rotochemical heating (see the upper panel of figure 1), these reactions are 
almost completely blocked due to the energy gap when rinpi < A. However, several 
reactions are opened when rfnpi > A. This is the most important effect of Cooper 
pairing in rotochemical heating since the quasi-steady state will be reached when the 
restoring mechanism given by the neutrino reactions becomes sufficiently important to 
conterbalance the effect of the spin-down forcing mechanism. 

RESULTS AND DISCUSION 

Fig. 1 shows that, for the superfluid case, the reactions are blocked until the chemical 
imbalances overcome the value of the energy gap of the neutrons A„, as argued above. 
These chemical imbalances are higher than those achieved in the non-superfluid case. 
This effect lengthens the timescale at which the system reaches the quasi-steady state, 
and implies that, in the presence of superfluidity, the chemical energy is larger and dis- 
sipated later to reheat the star. This makes it possible to fit the observation of the MSP 
J0437-4715 [2], unlike the non-superfluid case. Finally, for several EOS in which mod- 
ified Urea reactions are the dominant neutrino emission mechanism, this observation 
constrains the energy gaps to lie in the range 0.05[MeV] < min (A„ + 3Ap, 3A„ + A^) < 

0. 45 [MeV] [5]. 
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